Tunnelling in quantum superlattices with variable lacunarity 
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Quantum fractal superlattices are microelectronic devices consisting of a series of thin layers of two 
semiconductor materials deposited alternately on each other over a substrate following the rules of 
construction of a fractal set, here, a symmetrical polyadic Cantor fractal. The scattering properties 
of electrons in these superlattices may be modeled by using that of quantum particles in piecewise 
constant potential wells. The twist plots representing the reflection coefficient as function of the 
lacunarity parameter show the appearance of black curves with perfectly transparent tunnelling 
which may be classified as vertical, arc, and striation nulls. Approximate analytical formulae for 
these refiection-less curves are derived using the transfer matrix method. Comparison with the 
numerical results show their good accuracy. 



I. INTRODUCTION 

Quantum superlattices are microelectronic devices 
composed of alternating layers of two semiconductors 
with different energy band gaps deposited over a sub- 
strate [H, 0, H, The sandwiching of two materials 
whose band gap energies are different results in a re- 
tangular potential well for the electrons (holes) in the 
conduction (valence) band. The thickness of each semi- 
conductor layer is about 10-100 atoms (1-10 nm), so this 
three-dimensional finite potential well may be split into 
a two-dimensional system in the transversal dimensions, 
where electrons behave as a two-dimensional electron gas 
with continuous energy states, and a one-dimensional fi- 
nite square well in the growth direction with discrete 
bound states. 

In the absence of external electric field, semiconduc- 
tor superlattices may be considered, in the first approx- 
imation, as quasi-one-dimensional system of rectangular 
quantum wells separated by potential barriers (in the 
conduction band). When a external force is applied to 
the electrons in the quantum-well heterostructure they 
do not respond as free particles due to the influence of 
the other atoms. These effects may be taken into account 
by introducing an "effective mass" (m* ) which relates the 
particle motion to an external force (or potential) with- 
out worrying about all the atomic forces. Therefore, the 
scattering in semiconductor superlattices can be easily 
solved by means of the transfer matrix method using an 
effective mass for the electrons [1, H, H B • 

Fractal superlattices are designed by alternating semi- 
conductor materials according to an iterative fractal 
process using techniques such as molecular beam epi- 
taxy p^ . These devices may be crystalline or amor- 
phous. A typical device for the first kind is made by sev- 
eral layers of gallium arsenide (GaAs) of a few Angstrom 
(A) thick sandwiched between aluminium gallium ar- 
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senide (Ala:Gai_a;As) ones, where the crystal composition 
X may vary between and 1, such as those fabricated by 
Axel and Terauchi [TT| using X — 0.3, corresponding to a 
monolayer of Alj;Gai„2;As with thickness of 0.28 nm. A 
typical amorphous device is made by alternating layers of 
amorphous germanium (a-Ge) and amorphous silicon (a- 
Si) deposited in a silicon substrate [l3| , each layer having 
a thickness of less than 1.4 nm with good sharpness. 

The scattering of electrons in symmetrical polyadic 
Cantor fractal potentials, characterized by a lacunarity 
parameter which is independent of their fractal dimen- 
sion, have been numerically obtained by means of the 
transfer matrix method in Ref. [l^l- The representation 
of the reflection coefficient as a function of the particle 
energy and the lacunarity shows perfectly transparent 
tunneling along some curves, referred to as nulls, which 
may be classified as vertical (lacunarity- independent), 
arc (concave upward) , and striation (concave downward) 
nulls Ti^ . Analytical expressions for the position of these 
nulls have been obtained in Ref. |lj| in an optical context, 
which may be properly adapted to the present applica- 
tion. However, those authors omit the details of their 
derivation which may be of great interest when quan- 
tum fractal superlattices are incorporated as laboratory 
assignments [l^, ■ 

In this paper, the transfer matrix method is used to ob- 
tain analytical expressions for the reflection-less energies 
of particles in fractal superlattices designed as polyadic 
Cantor sets with variable lacunarity. The contents of this 
paper are as follows. The next section recalls the main 
facts about these kind of fractal superlatices. Section 3, 
presents the analytical formulas for the reflection-less en- 
ergies, whose validity is shown in Section 4 by its com- 
parison with numerical results. Finally, the last section 
is devoted to the main conclusions. 



II. POLYADIC CANTOR SUPERLATTICES 

In Cantor superlattices the distribution of the layers' 
thickness along growth axis corresponds to a flnite or- 
der approximation to a polyadic, or generalized. Cantor 
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FIG. 1: First two stages {S — 1 and S' = 2) of polyadic 
Cantor fractal sets witli A'^ = 5 (a) and A'^ = 6 (b), sliowing 
tlie definition of both the scale factor (7) and the lacunarity 
parameter (e). Black and white regions denote the potential 
values —V and 0, respectively. 



set [16|, |l7|. These sets are defined as follows. The first 
step [S = 0) is to take a segment of unit length, referred 
to as initiator. In the next step, 5* = 1, the segment is re- 
placed by N non-overlapping copies of the initiator, each 
one scaled by a factor 7 < 1. For even TV, as shown in 
Fig. [Tfb), one half of the copies are placed completely to 
the left of the interval and the other half completely to 
its right, each copy being separated by a fixed distance, 
let us say e. For odd iV, as shown in Fig. [U^a), one copy 
lies exactly centered in the interval, and the rest ones are 
distributed as for even A'', i.e., \_N/2\ copies are placed 
completely to the left of the interval and the other [iV/2j 
copies wholly to its right, where [TV/ 2 J is the greatest 
integer less than or equal to N /2. At the following con- 
struction stages of the the fractal set, S — 2, 3, . . ., the 
generation process is repeated over and over again for 
each segment in the previous stage. Strictly speaking 
the Cantor set is the limit of this procedure for S' ^ 00, 
which is composed of geometric points distributed such 
that each point lies arbitrarily close of other points of the 
set, being the S'-th stage Cantor set usually referred to 
as a pre-fractal or physical fractal. 

Symmetrical polyadic Cantor fractals are characterized 
by three parameters, the number of self-similar copies 
TV, the scaling factor 7, and the width of the outermost 
gap at the first stage, e, here on referred to as the la- 
cunarity parameter, as in most of the previous technical 
papers dealing wi th p olyadic Cantor fractals in engineer- 
ing applications [ij, ^M,- The similarity dimension of 
all polyadic Cantor fractals is D = In(Ar)/ ln(7^^) = 
— In(iV)/ ln(7), which is independent of the lacunarity 
parameter. The three parameters of a polyadic Cantor 
set must satisfy certain constraints in order to avoid over- 
lapping between the copies of the initiator. On the one 
hand, the maximum value of the scaling factor depends 
on the value of N , such that < 7 < 7max = 1/TV. On 
the other hand, for each value of N and 7, there are two 
extreme values for e. The first one is Emin = 0, for which 
the highest lacunar fractal is obtained, i.e., that with the 
largest possible gap. For even N , the central gap has a 
width of 1 — iV 7, and for odd iV, both large gaps sur- 
rounding the central well have a width of (1 — Nj)/2. 



The other extreme value is 

l-iV7 
N-2 ' 

1 -iV7 
N-3 ' 



even N, 



odd TV, 



where for even (odd) N two (three) wells are joined to- 
gether in the center, without any gap in the central re- 
gion. The width of the N — 2 gaps in this case is equal 
to Emax- Thus the corresponding lacunarity is lower than 
that for e = 0, but not the smallest one, which is ob- 
tained for the most regular distribution, where the gaps 
and wells have the same width at the first stage {S = 1) 
given by 

_ 1 - iV7 

Note that < Eiog < Cmax- 

Polyadic Cantor superlattices have been experimen- 
tally grown by Jarrendahl et al. (l^ by means of prop- 
erly alternating layers of a-Ge and a-Si. X-ray diffraction 
techniques were used to analyze the scattering problem 
of electrons in structures such a S = 5 triadic Cantor pre- 
fractal having 3^ (729) "elementary" alternating layers. 
A theoretical analysis of the diffraction peaks has shown 
good agreement with the results of the experiment. How- 
ever, a complete analysis requires taking into account 
both sample imperfections and instrumental limitations, 
which can influence location, intensity, and profile of the 
peaks. In this paper, in order to simplify the analysis, 
these effects are neglected. 



III. ESTIMATION OF THE NULL POSITIONS 

The graphical representation of the scattering reflec- 
tion coefficient as a function of both the particle energy 
and the lacunarity parameter for polyadic Cantor pre- 
fractal potentials, referred to as a twist plot, presents 
a characteristic null structure (black lines and curves) 
corresponding to energies at which the particle transpar- 
ently tunnels through the fractal potential wells, i.e., for 
which the reflection coefficient is null. 

The reflection coefficient nulls may be separated in 
three main families. The first and most noticeable is 
the family of vertical nulls, being caused by the bound 
states of the particle in each potential well, being thus 
independent of the lacunarity parameter. The nulls in 
the second family, denoted as arc nulls, are concave up- 
ward and curve from the upper left of the twist plots to 
the lower right. The third family of nulls, referred to as 
striation nulls, form the finer structure of the twist plots 
being curves which are concave downward, running from 
the lower left to the upper right of each twist plot. 

The nulls may be calculated from first principles as pre- 
viously done in Ref. [14] for the interference of light in 
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dielectric fractal superlattices, where only the final math- 
ematical expression is presented, completely omitting the 
details of the derivation. In the next subsections this gap 
is filled for quantum fractal superlattices. 



A. Transfer matrix method 

The scattering problem for polyadic Cantor sets with 
variable lacunarity can be easily solved by means of the 
transfer matrix method using an effective mass for the 
electrons. Let us summarize the main details of this 
widely known technique. 

The quantum scattering of particles on one- 
dimensional potential wells is governed by the steady- 
state, linear Schrodinger equation 



2 ni* dx^ 



V{x)'4}{x) =EiP{x), 



(1) 



where ipix), m* and E are the particle wavefunction, 
effective mass and energy, respectively, h is Planck's 
constant, and V{x) describes a distribution of potential 
wells, from here on, taken as a piecewise constant func- 
tion. The wavefunction i/ji on the region where the po- 
tential constant value is Vi, is the addition of two plane 
waves, ipi{x) = ipt(x) + ij;~{x), given by 



^± g±i ki X 



where i — ^/—l, ki is the local particle momentum, and 
Af are integration constants to be determined by apply- 
ing the standard boundary conditions at the interfaces 
between successive wells. In this paper, a distribution of 
potential wells of the same depth is considered, so the 
potential Vi = ~V in the wells and Vi = outside them. 

The solution of Eq. ^ for a distribution of N constant 
potential wells is easily obtained by means of the transfer 
matrix method [1, 0, S, 0, , yielding 



A 



N+l 



(2) 



where 



N 



M = Yl A n{di) D-^ Dn+1, 



D,, = 



1 1 



P^id^) 








where di is the width of the i-th potential well. 

Both the reflection and transmission coefficients of the 
scattering of a quantum particle, incoming from the left, 
with the TV-well potential is determined by the coeffi- 
cients of the matrix M, 



An 



Mil Mi2 
M21 M22 



A 



N+l 





where no backward particle can be found on the right 
side of the potential, so A'j^r 1 1 = 0. Both the refiection 
and transmission coefficients |19j are given by 



R 



T ■ 



Un 



\M: 



211 



1 4+ I 



I Mi: 



|2 ' 



(3) 



respectively, since k^+i — kg. 



B. Vertical nulls 



The vertical nulls shown in the twist plots are the re- 
sult of bound states of the particle due to its resonance 
on every potential well of the Cantor superlattice. Since 
there are potential wells at the S'-th stage pre-fractal, 
the vertical nulls are very noticeable in the twist plot. Us- 
ing the transfer matrix method is easy to calculate their 
exact positions, as shown in the appendix of Ref. [l3l |. 
This derivation is repeated here for completeness of the 
presentation. 

Let us consider a quantum well at the S'-th stage pre- 
fractal, whose width is a = L7^, where L is the length 
of the initiator, separated from the next one by a dis- 
tance b — Le^ — a{e/"f)^. For future convenience, let 
us take a symmetrical configuration in which the quan- 
tum well is surrounded by two regions of width 6/2, 
which corresponds to Vb = V2 = 0, Vi = —V, di = a, 
and do = d2 = 6/2, in the notation introduced in Sec- 
tion [TTTXl The particle momentum at the three constant 
potential regions may be normalized as 



ko = k2 



V2m*E 

h 



(4) 



^2m* {E + V) _ 



where (f> and (j)y are the (non-dimensional) particle en- 
ergy and depth of the potential well, respectively. The 
total scattering matrix is given by 



M = Po(b/2) Di Pi (a) D^^ Do P2(b/2) 



Mil M12 
M21 M22 



which may be easily calculated to obtain 



Mil = M2*2 ^e''''"' { cos(aki) + i M_t^ sin(aA:i) 

2 fco fci , 



(5) 
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Mi2 = Af2i = 



ftp 



2 fco /ci 



sin(a fci), 



where the asterisk indicates complex conjugate. 

The vertical nulls correspond to bound states of the 
particle inside the potential well, being characterized by 
a null reflection coefficient, for which Eq. ([3]) yields 



sin(afci) = sin(y02 + (^) = 0, 
having countably infinite solutions given by 

z==l,2,... 



(6) 



Note that the position of the vertical nulls is indepen- 
dent of both the initiator length and the pre-fractal stage 
thanks to the normalization of the energy 



C. Arc nulls 

The arc nulls are caused by collective bound states of 
the particle in any of the 2 A^'^~^ copies of the series of 
[-/V/2J potential wells wholly on either the right or left 
side of each of the TV"^"^ copies of the initiator at the S'-th 
stage of the pre-fractal superlattice. Since these collective 
bound states depend on the distance between individual 
wells, they are a function of the lacunarity parameter. 

The Cayley-Hamilton theorem ^2^ states that every 
matrix M satisfies the equation given by its characteristic 
polynomial which for a 2 x 2 matrix is given by pm{^) = 
A^-tr (M) A-hdet(M), where tr (M) and det(M) are the 
trace (sum of the diagonal elements) and determinant, 
respectively, of the matrix. For the scattering matrix 
of a potential well, Eq. ([5]), which is an unitary matrix, 



det(M) = |Mii| 



Ml 



= 1 and tr (M) e R. Without 



loss of generality, tr (M) = 2 cos(6'), where is the Bloch 
phase, being a real number if |tr {M)\ < 2 and a complex 
one otherwise, being pure imaginary only if tr (M) > 0. 
Therefore, the Cayley-Hamilton theorem yields 

= 2 cos(6') M - /, 

where / is the identity matrix. This equation allows the 
calculation of the n-th power of the matrix M as 

^^„^ sinM) ^^_ sin((n- l)g) ^ 



sin(6i) 



sin(6') 



which can be easily proved by the induction principle, in 
fact, note that 

M»+i =MM- = (2 cos(^?)) M - 

sin(6') ^ ^ " sin(6l) 

sin((n - 1) 6*) 



sin(6') 



■M, 



and use the formula 2 cos(6') sin(ri6') = sin((n -I- 1) 6*) -I- 
sin((n-l) Q). The off-diagonal element (M")2i of Eq. 



IS 

sm^fl) 

hence the reflection coefficient is null for 

sin(n6l)=0, sin(6i) ^ 0, (8) 

whose countably infinite solutions are the set of Bloch 
phases given by 

nOij = {ni + j)7r, i = 0,1,2,..., j = 1, 2, . . . , n - 1, 

(9) 

where n = [iV/2j in our case. Note that, in the right- 
hand side of this expression, tt appears multiplied by all 
the natural numbers except the multiples of n, for which 
sin(n6')/sin(6') — n. 

The Bloch phase may be calculated from Eq. ^ which 
yields 



tr (Af) — 2 cos(6fco) cos(afci) 
^0 ~^ ^1 



hi ki 



sin(6fco) sin(aA:i) 



2 fco fci 



cos(6 fco + a fci) 



_ (fco- fci )^ 
2 fco fci 

Equation (|4]) yields ak^ — <j) and 



cos(& fco — a fci). 



a fci = y 0^ + = 

so for (j) 3> (pv , fco ~ fci and 
tr (Af) = 2 cos(6l) 




2 cos(afci + 5fco) + 0((fci - fco)^) 



(10) 



Therefore, from Eq. ([9]), the normalized energy (pij for 
the arc nulls may be calculated as function of the lacu- 
narity as 



] 



z = 0,l,2,..., j = l,2,...,LiV/2j-l. (11) 
Note also that 

s\n{n6) = sin(n arccos(tr (Af)/2)) 



v/4-tr(Af)2 



C/„(tr (A^)/2), 



where Un{x) is the n-th Chebyshcv polynomial of the 
second kind 7]. The null value of this expression at the 
arc nulls is straightforwardly checked since 

tr(Af) = 2 cos(6',y) 2 (-1)^ 

where Eq. ^ has been used. 
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D. Striation nulls 

The striation nulls are caused by the quantum interfer- 
ence between the \_N/2\ potential wells on the right with 
those on the left side in each of the N^~^ copies of the 
initiator at the 5'-th stage of the pre-fractal superlattice. 
This interaction is exact for even N and only approxi- 
mate for odd N due to the presence of the potential well 
at the center of each copy of the initiator. Obviously, 
they also are a function of the lacunarity parameter. 

Let us start with the even N case. Let c be the distance 
between the rightmost potential well of the [-/V/2J wells 
on the left and the leftmost one of those at the right 
in every copy of the initiator at the S'-th stage, yielding 
c — L {1 — N — [N —2) e^). The scattering matrix for 
this problem is where 

M = Po(c/2) Af L^/2J Po(c/2). 

Equations ([5]) and ^ ior n — 2 yield 

2^, = (2i + l)7r, i = 0,1,2,..., 

where tr {M) — 2 cos{d), which may be easily calculated 
from 

tr (M) =tr (ML^/2J)tr (Po(c)). 
Introducing the trace operator into Eq. ([7]) yields 

sin((n - 1)61) 



sm 

tr (M") = . \ ' tr (M) - 2 



= 2 cos(n6'). 



sin(6i) ' ' sm{e) 
hence, using Eq. (fTO| . for cj)^ (jjy^ 

tr (M) = 2 cos(n {aki+bka)) 2 cos(cA;o)+0((fci - fco)^) 

yielding 6 k, n a ki + [nb + c) , where n = L-^/2j , 
resulting in 



\_N/2\ + + 



+ U-s _ TV - (TV _ LAr/2j - 1) 



(2z + l)-, z = 0,l,2,. 



(12) 



for the striation nulls for even N . 

For the odd N case, the central potential well between 
both series of [-/V/2J wells makes their mutual interfer- 
ence not exact. However, an approximate formula may be 
obtained as easily as before. Let c — L{1 — Nj^ — N e^) 
such that c/2 -I- 6 be the distance between the central po- 
tential well and that in the right (or left) in every copy 
of the initiator at the S'-th stage. The corresponding 
scattering matrix is Af ^ where 

M = Pi (a/2) D^^ Da Pa{c/2) A^ L^/^J • 
■Po{c/2)D^'D,P,{a/2), 



where the central potential well has been divided in two 
parts, one located completely to the left and the other one 
completely to the right. By exactly the same arguments 
as before, omitted here for brevity, tr (Af) ~ 2 cos{9), 
where 



{N~ lN/2\) (aki+bko) +cko, 



resulting in 



{N-lN/2\) ^<t>l+^l + 
+ (7~^-iV-(LiV/2j-l) 



= {2i- 



i = 0,1,2,, 



(13) 



as an approximation for the striation nulls for odd N . 

Note that Eg. (f T2 |) is a particular case of Eq. (fT3|l . 
hence in Ref. [ij] only the last one is shown for both 
even and odd N. 



IV. PRESENTATION OF RESULTS 

Twist plots are the gray-scale representation of 
the scattering reflection coefficient (in decibels, i.e., 
10 logj^Q R) as a function of the normalized energy (0) 
and the lacunarity parameter (e). Here, a linear gray 
scale was used, from black for null values to white for the 
maximum value equal to unity. Figures [2] (a) and [3] (a) 
show twist plots for, respectively, pentadic [N — 5) 
and hexadic {N = 6) Cantor fractal potential both 
with S = 1, 7 = 1/7, and <j)v = 1/2. The horizon- 
tal interval covers the normalized energy (p from to 

(3 7r)2 — 0y « 9.41, the energy of the third vertical 
null line. These figures have been obtained by a straight- 
forward numerical implementation of the transfer matrix 
method. 

The most noticeable fact in Figs. [5] (a) and [3] (a) is 
the null structure (black lines and curves) corresponding 
to energies at which the particle transparently tunnels 
through the fractal potential wells. The strongest (dark- 
est) among these curves are the vertical nulls. The arc 
nulls, which are concave upward, curving from the up- 
per left to the lower right of the plots, appear in pairs 
for iV = 6 (see Fig. [3] (a)), as expected from Eq. (fTTj) . 
The striation nulls, which are concave downward, run- 
ning from the lower left to the upper right of each twist 
plot, form the finer structure of the twist plots. 

The main difference between Figs. [D (a) and [3] (a) is 
observed for iV = 5 and, in general, for odd N , because 
in such a case the tunnelling is not perfectly transparent 
at the whole length of the arc nulls, but instead appear 
dashed, alternatively darker and lighter after they cross 
the striation nulls. The striation nulls also appear lighter 
when they cross the arc nulls, resulting in a leopard- 
like spot pigmentation pattern. This result is not un- 
expected, since both the arc and the striation nulls are 
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FIG. 2: Scattering reflection coefficient (a) for tlie poiyadic 
Cantor pre-fractal potentiais witli S = 1, N — 5, y = 1/7, and 
(j)v = 1/2, and tfie analytical curves (b) for the correspond- 
ing vertical (continuous line), arc (dashed line), and striation 
(dotted line) nulls. 



FIG. 3: Scattering reflection coefficient (a) for the poiyadic 
Cantor pre-fractal potentials with S = 1, A*' = 6, 7 = 1/7, and 
(l)v = 1/2, and the analytical curves (b) for the correspond- 
ing vertical (continuous line), arc (dashed line), and striation 
(dotted line) nulls. 



V. CONCLUSIONS 



not truly null for odd N due to the single quantum well 
between each pair of lN/2\ potential wells, which intro- 
duces non-periodicity in the scattering problem resulting 
in a (small) non-null reflection coefficient. 

The nulls calculated by Eqs. (HI]), and ((13]), are 
shown in Figs. [D (b) and [3| (b). These figures clearly 
show the good accuracy of the analytical approximations 
which reproduce the main features of the corresponding 
twist plots. For even N the accuracy is surprisingly good. 
For odd N, as expected, both arc and striation nulls are 
less accurate. Note also that the accuracy degrades for 
small energy as expected due to the lack of validity of the 
approximations used in Sections IIII CI and If 11 Dl in such 
a case. 



Analytical expressions for the calculation of the parti- 
cle energy for transparent tunnelling in poiyadic Cantor 
pre-fractal potentials as a function of the lacunarity pa- 
rameter have been obtained from first principles by using 
the transfer matrix method. The comparison with the 
results obtained by a numerical implementation of the 
transfer matrix method shows the good accuracy of these 
expressions for A^-adic Cantor pre-fractals with even N 
and their reasonable accuracy for odd N. The reasons 
for this difference in accuracy have also been presented 
in detail. 

The analytical results obtained in this paper can be 
easily incorporated as a complement in computer labo- 
ratories for undergraduate quantum mechanics or solid- 
state physics courses currently using the transfer matrix 
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method only as a numerical tool [H, [TB|- In fact, the 
possibility of approximately calculating the geometrically 
complex figures observed in the results provides the stu- 
dent with an excellent example of the combination of 
theory and numerical experiments. Thus, different as- 
pects of the model may be assigned to different students 
or groups of students. In fact, stimulating discussions 
among the students follow naturally from this kind of 
approach. 

Finally, since the transfer matrix methods may be 
adapted by physics analogy to the propagation of waves 
in general one-dimensional quasiperiodic media [2]| . the 
present results may be also applied to, for example, 
acoustics, optics, or vibrating strings. In fact, the effect 
of the lacunarity in the scattering on generalized Cantor 
media was first studied in optics [H, [3| , where the ad- 



ditional degree of freedom introduced by the lacunarity 
was used with success to obtain new physical character- 
istics. A detailed study of these applications, outside the 
scope of this paper, may also be an estimulating topic for 
further study. 
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